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Abstract. All complete, axially symmetric surfaces of constant mean cur- 
vature in R 3 lie in the one-parameter family D T of Delaunay surfaces. The 
elements of this family which are embedded are called unduloids; all other 
elements, which correspond to parameter value r £ R — , are immersed and 
are called nodoids. The unduloids are stable in the sense that the only global 
constant mean curvature deformations of them are to other elements of this 
Delaunay family. We prove here that this same property is true for nodoids 
only when r is sufficiently close to zero (this corresponds to these surfaces 
having small 'necksizes'). On the other hand, we show that as r decreases to 
— oo, infinitely many new families of complete, cylindrically bounded constant 
mean curvature surfaces bifurcate from this Delaunay family. The surfaces in 
these branches have only a discrete symmetry group. 



1. Introduction 

In 1841, C. Delaunay discovered a beautiful one-parameter family of complete 
noncompact surfaces of constant mean curvature one in K 3 which are invariant 
under rotations about an axis @. (Henceforth we shall abbreviate constant mean 
curvature one by CMC.) Using the rotational symmetry the search for these sur- 
faces reduces to finding their meridian curves, and these may be found, in turn, by 
solving an appropriate ODE. This family of 'Delaunay surfaces' is parametrized by 
a variable r lying in the set 

V= (-oo,l] -{0} =V + UV-; V- = (-oo,0), V + = (0,1]. 

The Delaunay surface corresponding to any value r £P will be denoted D T . When 
r G V + , then D T is embedded and is called an unduloid; the meridian curve in this 
case is the roulette of an ellipse. When r£P^, then D T is no longer embedded, and 
is called a nodoid; the meridian curve of any nodoid is the roulette of a hyperbola. 
There is a good geometric limit of the surfaces D T as r — * (which is the same 
for r \ and r / 0), which is the nodded surface formed by the infinite union 
of mutually tangent spheres of radius 1 arranged along a common axis. A nice 
geometric description of these Delaunay surfaces can be found in [H . 
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In this paper we shall prove the existence of some rather surprising new families 
of complete CMC surfaces which arise as deformations of Delaunay nodoids. Let us 
digress briefly to set these new surfaces into their proper context before describing 
our results in more detail. 

The past several years have witnessed great progress in the understanding of 
more general complete, CMC surfaces of finite topology. On the one hand, the 
great variety and flexibility of such surfaces has been suggested by numerous nu- 
merical and computer experiments, notably by the GANG group at the University 
of Massachusetts and Grosse-Brauckmann and his collaborators at the Universitat 
Bonn, see jittp : / / www. gang . umass . edu / cmc] for many nice pictures. On the other 
hand, there have been several theoretical advances which vindicate many of these 
computer experiments. These advances may be subdivided into two groups of re- 
sults, those concerning various constructions of complete CMC surfaces, which are 
not relevant to the present discussion, and those concerning the structure theory of 
these surfaces. 

This structure theory is much better developed when we restrict attention to 
complete (oriented) Alexandrov-embedded CMC surfaces of finite topology, the key 
point being that any one of these surfaces has only finitely many ends, and each 
of these ends is necessarily modeled on a Delaunay unduloid. Recall first that 
any surface of this type is conformally equivalent to a punctured Riemann surface 
£ = £ \ {pi, . . . ,pk}; let Y 3 be a handlebody such that dY = £. Then £ is said 
to be Alexandrov embedded if its immersion into R 3 extends to an immersion of 
Y \ {pi, . . . ,pk}- In particular, D T is not Alexandrov-embedded when r < 0. 

A deep theorem of Meeks 1 14 ] states that any end of an Alexandrov-embedded 
CMC surface of finite topology is cylindrically bounded. Following this, Korevaar, 
Kusner and Solomon [Q showed that for any such end there is a Delaunay unduloid 
D T to which the end converges exponentially. contains a further global struc- 
ture theorems. These results show that Delaunay unduloids play a fundamental 
role as building blocks for finite topology Alexandrov-embedded CMC surfaces. In 
addition, this tameness of the geometry of the ends leads to a fairly detailed un- 



derstanding of many analytic and geometric problems on these surfaces, M, Jl2|. 

It is natural to try to extend this theory to complete, CMC surfaces with finite 
topology which are not necessarily Alexandrov-embedded, for example those with 
ends modeled on Delaunay nodoids. Fairly general examples of such surfaces are 



constructed in |11|, [12], but the ends of those surfaces are modeled on Delaunay 
nodoids with r very close to zero. 

The results of this paper begin to elucidate the limitations of the structure 
theory Q, and how it fails when Alexandrov-embeddedness is dropped. More 
specifically, we examine the stability theory of the Delaunay nodoids D T as r — > 
— oo. In a sense we explain below, there is some stable range r» < r < for 
which the only deformations of D T are to other Delaunay surfaces. However, as 
r decreases, we establish the existence of new CMC surfaces which bifurcate from 
the Delaunay family. These surfaces have only a discrete symmetry group, rather 
than full rotational symmetry, and the number and geometric complexity of these 
bifurcations increases as r becomes more negative. All of these new surfaces remain 
cylindrically bounded, but their geometric structure is rather intricate. 

To conclude this general discussion, we recall that the basic theory of the moduli 
space of complete Alexandrov-embedded CMC surfaces of finite topology is estab- 
lished in 0. The paper [12| extends this theory to include CMC surfaces with ends 
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modeled on Delaunay nodoids with r > r* ; (it also develops some global aspects of 
this moduli space theory) . Our results here show that a comprehensive description 
of the moduli space theory of complete immersed CMC surfaces of finite topology 
must incorporate these new families of CMC surfaces. 

We may now give the precise statements of our main results. These require a 
few preliminary definitions. 

Definition 1.1. We set terminology for a family of rigid motions in R 3 , nor- 
malized for convenience to fix the z-axis, as well as the surfaces which have the 
associated symmetries. 

i) Let j £ N, j > 2. Let Rj denote the rotation by angle 2ir/j about the z-axis 
in R 3 . We say that a surface E C R 3 is R } ; -symmetric if it is invariant by 
Rj, but not by any Rji for j' > j. 

ii) For a £ [— tt, tt], let R a denote the rotation by angle a in R 2 . Also, let 
t a £ R. We say that the rigid motion S a of R 3 is a screw motion about the 
z-axis (of angle a and translation length t a ) if it has the form 

S a : R 2 x R — > R 2 x R 

(xi,x 2 ,z) i — > (Ra x, z + t a ) 

We say that the surface S C R 3 is S a -symmetric if it is invariant with respect 
to some screw motion S a , for some choice of translation length t a , and if it 
is not invariant with respect to any other screw motion S a i with < a' < a 
if a £ (0, tt] or a < a' < if a £ [-tt, 0). 

iii) Finally, we say that the surface S C R 3 is Tj. Q -symmetric if it is both Rj- 
symmctric and S Q -symmetric. 

Our first main result states that there are infinitely many bifurcations from 
the Delaunay family D T as r — > — oo to new CMC surfaces which have a discrete 
symmetry group: 

Theorem 1.2. Let j > 2; then therefor any a £ [— ir/j, it/j] — {0}, there exists 
a Tj^ a - symmetric CMC surface. In fact, each of these surfaces lies in a branch which 
bifurcates from the nodoid D Tj at some value Tj < — \J p — 2. 

This may be proved either using the general bifurcation theorem of Smoller and 



Wasserman [17| or more simply by degree theory [16|. As in most such problems, 
the main point is to show that the index (i.e. number of negative eigenvalues) of 
the Jacobi operator L T on some compact quotient of D T increases as r \ — oo. As 
part of this, we show that there exists some value t» £ (— 2, — v2) such that no 
bifurcations occur when t* < r < 0. 

This result can be sharpened when j is large enough, for in this case we can 
guarantee that the bifurcation value Tj is large negative. This allows us to obtain 
better control on the spectrum of Jacobi operator, and so we may use the bifurcation 
theorem of Crandall and Rabinowitz |2j to show that the bifurcation branches are 
smooth. 

Theorem 1.3. There exists jo > 2 such that when j > jo and a £ [— Tt/j, Tt/j] — 
{0}, then there are values Tj^ a < T)j }0! > 0, and a real analytic branch of Tj, a ~ 
symmetric CMC surfaces, which we denote Dj. a (T,rj), parametrized by r\ with |ry| < 
nj^ a , which bifurcate from D T .. These families of surfaces have the following 
properties: 
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(i) As 77 — s- 0, the surfaces Dj ta {rf) converge uniformly on compact subsets to 
the nodoid D T . a , in C°° topology. 

(ii) Locally on D T . and for n sufficiently small, we can write Dj t0l {rj) as a 
normal graph over D Tj a for some function . The functions w^/n converge 
uniformly on compacts sets to a nontrivial function of the form (s, 9) — > 
4>(s) cos(j9). (The variables (s,9) will be defined in $2.) Furthermore, the 
function ip(s, 9) := 4>(s) e %e is a nontrivial solution of C T ip — which satisfies 
cf>(s + 2irs T ,9) = e ila 4){s,9). 



In fact, we prove that 



Remark 1.4. It is likely that Theorem 1.3 holds whenever j > 2. We also 
suspect that for any j > 2, \a\ < ir/j, a 7^ 0, there are precisely two bifurcating 
branches, and these correspond to two separate bands of continuous spectrum of 
£ T , with one branch bifurcating at the value r = Tj, Q , which lies in the first band, 
and the other bifurcating at a different value Tj, T lying in the second band. This 
second bifurcation point should satisfy 

fj a = — j = + o{r l ). 

v/l-(|aj|/27r-l)2 

as j — > 00. 

We conclude this introduction by mentioning that the bifurcations we find are 
related to the so-called 'Rayleigh instability of the cylinder'. We refer in particular 
to jl], §111, for a discussion of this phenomenon as manifested in the capillary 
instability of a liquid jet. As an historical aside, this problem was originally studied 
by Plateau, who posited that a cylindrical jet should break up into rotationally 
symmetric pieces. Lord Rayleigh found fault with Plateau's argument but emended 
it by noting that in certain situations the instability should produce nonrotationally 
symmetric perturbations. 

Other types of 'CMC cylinders' with few or no symmetries are known to ex- 
ist. For example, it is pointed out in 0] that a 'relaxation' of the construction of 
Wente tori produces CMC surfaces with finite rotational and discrete translational 
symmetry, or with discrete screw motion symmetry. In addition, one can produce 
immersed CMC surfaces with continuous screw motion symmetry using the clas- 
sical associate family construction. Finally, [jlj] uses the DPW (Dorfmeister-Pedit- 
Wu) method (which is akin to the Weierstrass representation formula for minimal 
surfaces) to produce examples of cylindrically bounded CMC surfaces with no sym- 
metry. To our knowledge, none of these lie in a continuous family including the 
Delaunay surfaces. 

We are grateful to the referee for drawing to our attention the work of Lord 
Rayleigh, as well as the examples in the last paragraph. 

2. Isothermal parametrization of Delaunay surfaces 

Since the Delaunay surfaces D T are surfaces of revolution, their most natural 
parametrizations would seem to be the obvious ones: 

(u, 9) 1 — ► {p(u) cos 9, p(u) sin 9,u), 
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where function p — p T is a solution of an ODE which is derived from the constant 
mean curvature condition. However, for most analytic purposes it turns out to 
be far more convenient to use a different parametrization which is isothermal, and 
which we now describe. This definition may seem ad hoc, but is motivated by a 
systematic line of reasoning in the theory of integrable systems; [11] contains a 
detailed derivation of the fact that this parametrizes D T . 

This isothermal parametrization rests on two functions a = a T and k = n T , the 
definitions of which vary, according to whether r is positive or negative. First, a is 
the unique smooth nonconstant solution to the initial value problem 

(2.1) + 7-2 cosh2 <7 = 1 > ds(J ^ = °' < °' when T G V+ ' 

(2.2) f-^J +T 2 sinh 2 (T = 1, <9 s ct(0) = 0, <r(0) < 0, when r e V~ . 
Next, k(s) is the unique solution of 

(2.3) = r 2 e ff cosher, k(0) = 0, when r e "P+ , 
as 

(2.4) = -t 2 e CT sinha, k(0) = 0, when r E V~ . 
as 

The change of variables from the previous cylindrical coordinates (u, 9) is ef- 



fected by setting u — n(s) and one then has (see [11]) p(n(s)) — Te a ^ s \ and so the 
new isothermal parametrization is given by 

(2.5) X T :lx S 1 3 (s,9) i — > - (t e aAs) cos 9, r e° As) sinfl, k t (s)) . 

The metric coefficients in this new coordinate system are g ss = ggg = T 2 e 2<T , g s g = 
90s = 0. 

Regardless of the sign of t, the function s — > er T (s) necessarily changes sign. 
Together with (|2.4| ), this means that when r < 0, then k is not monotone; in 
contrast, using ( |2.3| ), k is monotone when r > 0. Hence although the differential of 
X T is always full rank, so that X T is always an immersion, it is only an embedding 



when t > 0. Notice also that when r = 1, then (2.1) implies that a = 0, and so 



k(s) = s and p = 1/2, which means that D\ is a cylinder of radius 1/2. 

3. The period function 



The Hamiltonian nature of the equations (2.1) and (2.2) implies that a is 
periodic. We denote its period by a 2tts t . We now investigate the dependence of 
s T on t. 

To begin, s T has an integral representation. To state it, let a T ,a' T > be 
determined by the equations r cosh a T = 1 when t > and r sinh a' T = — 1 when 
t < 0. Then 

(3.1) s T :=- / - when r£p+, 

n Jo VI - t 2 cosh 2 cr 

and 
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We also use the following equivalent forms of these representations 
2 f bT dx 

(3.3) s T :=- / when reP+, 

Jo V cos z x — r z 

where fe T = arccosr € [0, 7r/2), < r < 1, and 
2 dx 

(3.4) s T :=- / - when reF. 

These are effected by the changes of variables r cosher = cosx and rsinh a — cosx, 
respectively. 

Proposition 3.1. As a junction ofr, s T is monotone increasing when — oo < 
t < and monotone decreasing when < r < 1. 

PROOF. It is obvious from (|j|) that ds T /dr > when t < 0. The fact that s T 
decreases when r increases from to 1 is less obvious and somewhat more difficult 
to obtain; since we do not need it here we omit the proof. □ 

We next consider the asymptotics of s T as r approaches various possible limiting 
values in V . The first two cases are less important, but we record them anyway; 
namely, we have 

(3.5) s T = --logr 2 +0(1) as r^O 



and also 



lim s T = 1. 



On the other hand, the behaviour of s T as r — * —oo is fundamental to our analysis. 
To study this we introduce an auxiliary function 

7 T (i) := toy (s T t) , t < 0; 
by design, 7 has period 2ir. 

Lemma 3.2. As r — > —00, 7 r converges uniformly to the function cost. More 
precisely, 

(3.6) j T (t) =cost + 0(\r\- 2 ) and d T j T (t) = 0(\t\- 3 ), 
as t — ► — 00 . Simultaneously, in the same limit 

(3.7) s t = -- + 0(\t\- 3 ), and d T s T = \ + G(|t|- 4 ). 



Remark 3.3. The proof below easily extends to give full asymptotic expan- 
sions for s T and 7 in powers of 1/r as r — ► —00. However, we only require the 
terms given in this statement. 

Proof. Set er := c(0). Since |ct| < — <r n , we may define the function w(s) by 

cr(s) = ct(0) cos(tw(s)). 



To normalize it, we require that w(0) = and w'(0) > 0. The equation for a, (2.2), 
becomes 

2 sinh 2 cr - sinh 2 cr 
H 5 = 0, 
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or equivalently 

(3.8) (dsw) 2 = - y ) where $(t) := sinh 2 VI. 



(2.2) also implies that \a\ < | sinh cr| < \t\ 1 . Using the first term of the Taylor 



expansion of $ gives d s w = 1 + C(|t| 2 ), and so 
(3.9) w( S )= S (l + 0(\rr 2 )). 

Next, the periodicity of a translates to the equality 

-2tt = tw(2its t ) = t(2tts t + 0{s t \t\- 2 )), 

or equivalently 

s t = --+0(\t\- 3 ). 

T 

Inserting this into the expression for 7(t), and using <jq — arcsinh(l/r) = 1/r + 
C(l/|r|- 3 ) gives 

7 (t) =cost + (D(\T\~ 2 ). 

This procedure may be continued, using this expansion for 7 (and hence for a) to 
determine the next term in the expansion for s T , then using this in turn to get the 
next term in the expansion for 7, and so on. □ 

4. Spectral analysis of the Jacobi operator 

The primary tool in this paper is a detailed analysis of the spectrum of the 
linearized mean curvature operator C T (which is usually called the Jacobi operator) 
on the family of Dclaunay surfaces D T . We begin by describing the specific form 
of this operator and then review how Bloch wave theory may be applied in this 
context. This theory reduces the spectral analysis to that of a countable collection 
of continuous families of operators with discrete spectrum and it shows that the 
spectrum of C T is a union of 'bands' of absolutely continuous spectrum. Geometric 
considerations allow us to identify specific solutions of C T u — 0, and this, in turn, 
allows us to track the location of some of these bands as r — > — 00. 

4.1. Coordinate expression. It is well known that the Jacobi operator has 
the form 

C T := — A T — |A T | 2 , 

where A r and A T are the Laplace operator and second fundamental form on D T , 



respectively. In terms of the isothermal parametrization (2.5), 
(4.1) C T = --^ {pi + d 2 e - r 2 cosh2a) . 

We shall be using the periodicity of the coefficients of this operator in a central 
way. However, it is more convenient to work with a family of operators for which 
the period is fixed as r varies. Accordingly we change variables, setting 

so that the coefficients of L T with respect to (t, 9) have period 2ir for all r S V . In 
these new coordinates, 

C r = - taiaj (df + 4d 2 e - s 2 T r 2 cosh(2 CT )) . 



8 



RAFE MAZZEO AND FRANK PACARD 



Removing the factor l/(s 2 T 2 e 2<T ), we define 

L T := -d 2 - s 2 T dg - s 2 t t 2 cosh(2cr). 

Since we are only interested in some aspects of the spectral analysis of C T , such 
as the number of negative eigenvalues, the existence of a nullspace, etc., we may 
concentrate on the study of the slightly simpler operator L T . 

4.2. Spectral decomposition. The operator L T has many symmetries, and 
these may be used to reduce it to operators for which the spectral analysis is more 
tractable. 

The first and most obvious reduction uses the rotational invariance in 6. Thus, 
if u(t, 9) £ L 2 (M. x S 1 ), we have the decomposition 

u(t,0) = ^ Ui (t)e yfl , 

iez 

where each of the coefficients Uj(t) is in L 2 (R) and X) ll u ill 2 < 00 • The operator 
L T induces the operator 

(4.2) L Ttj = -8 2 + s 2 T j 2 - s 2 T t 2 cosh(2<r), 

on the j th eigenspace. In fancier language, the operator L T is reduced by the 
splitting 



L 2 (Rx S 1 ) =0£ 



from this eigenspace decomposition into a direct sum of (self adjoint) ordinary 
differential operators L T j. Note that L T j — i T ,-j, even though the spaces L 2 and 
L 2 _j are different (one corresponds to the eigenfunction e^ 6 and the other to the 
eigenfunction e"^ 6 ). It is clear that 

spec(L r ) = I^J spec(L TJ ). 

Furthermore, noting that 

(4.3) L T:j = L Ti0 + s 2 T j 2 , 

we may as well restrict attention to the operator L T and the single Hilbert space 
L 2 (R). 

To analyze this last operator we use the technique of Bloch waves, also known as 



Flocquet theory. We refer to [15 1 and 10] for details, cf. also [13]. This relies on a 
direct integral decomposition of L 2 (R) defined using the Fourier-Laplace transform. 
Given any function / in the Schwartz space S on R, define 

jez 

We may as well assume that — it < a < it. This function satisfies f(t + 2tt, a) — 
e la f(t,a), and so we lose no information by restricting t to lie in [— 7T,7r]. It is 
straightforward to check that 

H/|li 2 (H) = 27r ll/llL 2 ([-^,7r]2), 
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which means that the map / i— > f may be extended as an isometry from L 2 (R) into 
L 2 ([— 7r, 7r] 2 ). In fact it is an isometry since there is an inversion formula given by 

f(t) = ^- r f(i,a)e ia Ua, 

where i e [— tt, it] and t = i + 27rj. Altogether, this gives the direct integral 
decomposition 

J Q£[- 7T,7r] 

the Hilbert spaces L 2 a appearing here are all the same, but the index reflects the 
fact that they arise as the closures of the space of continuous (or H 1 ) functions 
u which satisfy u(ir) = e la u(—Tr). The reason for this distinction will be evident 
soon. 

The importance of this decomposition for us is that the operator L Ti0 reduces 
further into a direct integral of operators i?o,a( T ) induced on each one of the spaces 
L 2 a . More specifically, Hq i01 (t) has the same expression as L Tt o, but its domain is the 
space of H 2 functions on [—w, n] satisfying the quasiperiodic boundary conditions 
u{tt) = e M u(-7r), u' {it) = e ia u'(-ir). 

As an aside about notation, we will be considering the reduction of the full 
operator L T to the subspace of a-quasiperiodic functions, and shall denote this op- 
erator by H a (r). Its restriction to the j th eigenspace of the cross-sectional operator 
is denoted Hj ta (r), and since 

(4-4) H 3 . a (T)=H , a (r)+4f, 

it will suffice at first to consider only the case j = 0. Notice that the subscripts 
j and a refer to the symmetries of this reduced subspace (i.e. j th eigenvalue on 
the cross-section and quasiperiodicity a); the parameter r has been elevated from 
subscript status to ground level, since our ultimate concern is with the spectral flow 
of these operators as a function of r. 

Continuing the main thread of the discussion, the spaces L 2 a are all the same, 
but the subscript a is meant to remind us that the domains of the restrictions of L Tj o 
to each one of them is different. Since these are self-adjoint boundary conditions, 
each H 0ya (T) has discrete spectrum: 

spcc(i7o,a(r)) = {A (r, a) < Xi(r,a) <...}. 

The k th eigenvalue map can be considered as a map: 

Afe(r, •) : S 1 — ► R; 

its image is an interval Bfe(r), which is called the fc th band; the union of all these 
intervals constitutes the spectrum of L Tt0 - 

spec(£ r , ) = (J B fe (T). 

feeN 

A bit more is known about these bands, and the 'band functions' Afc. First, since 
L t ,q is an ordinary differential operator, the space of solutions of (£ T ,o — A)</> = 
(with no growth restrictions) is precisely two-dimensional. Suppose that A is in 
the spectrum, i.e. A = Afe(r, a) for some k and a. The corresponding eigenfunction 
satisfies ip(t + 2n) = e ta Lp(t): since the coefficients of the operator are real, 4> 
is in the same eigenspace and satisfies <p(t + 2n) = e~ la cf){t). This implies that 
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Afe(r, —a) = Afe(r, a), and hence we may as well restrict < a < tt. In fact, it is 



known [10] that the band functions with even index are nondecreasing while the 



band functions with odd index are nonincreasing in this interval, so that 

A (r,0) < A (T,vr) < Ai(r,7r) < Ai(r,0) < A 2 (r,0)... . 

This means that the bands are nonoverlapping, except perhaps at their endpoints, 
and we have 

B 2k {r) = [A 2 fc(r,0),A 2fc (T,^)] 

-B 2 fe+i(r) = [A 2 fc+i(r,7r), A 2 fc+i(r,0)]. 

No band Bk(r) may reduce to a point for any r G V . The reason is simply 
that if this were the case, then the band function Afe(r, a) would be constant in a, 
say Afc(r, a) = A°(r), and hence there would be an infinite dimensional space of 
solutions to the equation L T $u = A^(r)w, which is at odds with the fact that L T $ 
is an ordinary differential operator. 

For any A G R, the space of solutions of L T ^u = Xu is two-dimensional. 
When A G UkenBkiT), then the eigenspace contains two quasi-periodic (and hence 
bounded) functions, whereas if A lies outside the spectrum, then the correspond- 
ing set of solutions has a basis {u±}, where u+ decays exponentially as t — ► — oo 
and grows exponentially as t — > oo, while u_ grows exponentially as t — > — oo and 
decays exponentially as t — > oo. 

The main task lying ahead is to understand the dependence on r of these 
bands B^(t). While it is certainly not possible to do this explicitly, we will use 
the special geometric solutions of L T u = 0, the existence and nature of which give 
some nontrivial information about the bands B^{t), h = 0, 1, 2. 

4.3. Geometric Jacobi fields. Solutions of C T u = (or equivalently, L T u = 
0) are called Jacobi fields; they correspond to variations of the Delaunay surface D T 
which preserve the mean curvature to second order. We are particularly interested 
in solutions which defined on all of D T which grow at most polynomially. In general, 
Jacobi fields may not correspond to actual CMC deformations, but as we describe in 
a moment, these correspond to global CMC deformations which arise from explicit 
geometric motions. 

The Bloch wave theory of the previous subsection, cf. [13 , (9| can be used to 



deduce the Fredholm properties and asymptotic expansions as t — > ±oo for solutions 
of L T u = 0. A consequence of this development is the fact that Jacobi fields have 
definite exponential rates of growth or decay. In particular, one proves that there 
is at most a finite dimensional family of Jacobi fields u which satisfy the bounds 
|u|e _e '*' < C for some (and in fact any) sufficiently small e. Somewhat remarkably, 
it is possible to identify every one of these temperate Jacobi fields because all of 
them correspond to explicit CMC deformations of D T . 

There are three types of global CMC deformations: those corresponding to 
translations of D T , which are a three-dimensional family, those corresponding to 
rotations of D T , which are a two-dimensional family, and a final one-dimensional 
family corresponding to changing the Delaunay parameter t. The infinitesimal 
variations corresponding to each of these span the six-dimensional space of temper- 
ate Jacobi fields. We now describe the distinguished three-dimensional subspace of 
these which are bounded. 

We temporarily revert to the isothermal (s, 9) coordinate system. 



BIFURCATING NODOIDS 



11 



• There is a Jacobi field 3>°(s) which corresponds to translating D T along its 
axis. It is obtained by projecting the constant Killing field (0, 0, 1) generating 
this translation over the normal vector field to D T . This yields 

:= d s a T . 

This Jacobi field is rotationally invariant and is also periodic in s, hence 
bounded. In addition, <J>° has precisely two nodal domains on the portion 
of D T where s <G [0, 2 tt s t ] . 

• The two Jacobi fields corresponding to translating D T orthogonally to its 
axis are given by projecting the constant Killing fields (1,0,0) and (0,1,0) 
on the unit normal field. These functions are again periodic in s; appropriate 
(complex) linear combinations of them have the form <f>^ 1 (s)e ± ^. By closer 
inspection, we actually see that 

^f 1 = cosher e ±ie when reV + , 

and 

Qf 1 = sinhtre ±ie when reP". 

Notice that cosh a does not change sign while sinh o has precisely two nodal 
domains when s e [0, 2 tt s t ] . 

Remark 4.1. There are three other geometric Jacobi fields, corresponding to 
rotating D T around its axis and varying the Delaunay parameter. These are not 
periodic but instead grow linearly, hence do not enter our considerations below. 
Therefore we omit further discussion of them. 

4.4. The spectral flow. We shall henceforth always reduce to a-quasiperiodic 
functions, and so it is no longer necessary to think in terms of the operator L T on 
the complete Delaunay surface D T , nor the band structure of its spectrum; rather, 
we focus on the operators H a (r), which are simpler to analyze because they have 
discrete spectrum. The ultimate goal now is to determine when these operators be- 
come unstable, i.e. when they have nonzero index (number of negative eigenvalues). 
We know this does not occur when r > 0, nor when < t < 0. In this section we 
analyze the spectral flow of H a {r) as t decreases to — oo. 

The eigenvalues of this operator take the form Xkjij) = Afe(r, a) + s 2 j 2 , j, k = 
0,1,2,..., and so this study divides into essentially two parts. First we obtain some 
information about the values of \k{r, a) when k = 0, 1, 2; following this we turn to 
the problem of when Afej(r) < 0. 

We take up the first of these tasks in the next two subsections. Because Afe(r, a) 
always lies between Afe(r, 0) and Afe(r, tt), we focus most of the attention on functions 
satisfying periodic or antiperiodic boundary conditions. 

4.4.1. The spectrum of Hq i(x {t) when r > 0. Although not needed later, we 
first show that there is no spectral flow for the operators H a (r) as r increases from 
to 1. 

Proposition 4.2. Suppose that t e V + . Then the bottom of first band B (t) 
occurs at A = — s 2 ., i.e. 

A (r,0) = -s 2 r . 

In addition, is either at the top of the second band or bottom of the third, i. e. 
either Ai(r,0) = or A 2 (t,0) = 0. 
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PROOF. We have indicated that L T (cosh ere 19 ) — 0, and it satisfies periodic 



boundary conditions, so that ffi j o(r)(cosh(7) = 0. By (4.4), this is the same as 

Hq,o(t) (cosher) = —s 2 cosher. 

But since cosher is everywhere positive, it must correspond to the ground state 
eigenvalue. 

Next, the function <3?° = d s a is a solution of L T w — 0, and again satisfies 
periodic boundary conditions; moreover, it has precisely two nodal regions in a 
period domain, and so it must correspond either to the second or third eigenvalue 
of H 0: o{ T ), i-e. either to Ai(r, 0) or A2(r, 0), as stated. Unfortunately, it is not clear 
from the evidence at hand which of these is actually the case. □ 



Corollary 4.3. When t > 0, the bands Bk{r), k > 3, are strictly contained 
in the positive half-line R + , while B^ir) C M + (possibly in the open half-line). In 
any event, the operator Ho t0t (r) has no spectral flow as t varies in V + . 

The first statement follows at once from the previous Proposition. Obviously, 
since Afe(r,a) > — s 2 ., we have \kjij) > for all j > 1, k > and all r € [0,7r]. 
This proves the last statement. 

4.4.2. The spectrum of Hq u (t) when r < 0. In contrast to the situation for 
unduloids, there is more to say about the spectrum of H a (r) for nodoids. 

Proposition 4.4. Let r e V~ . Then for a S [0,7r] we have 

- 2 »' + (s) 2 - A ^ A »<-°> s (S 2 - A ' ; 

in addition, 

A 1 (r,0) = -«», A 2 (r,0) = 0. 



Proof. First consider the identities for Ai and A 2 . We already know that 
L r (sinh<7 e l9 ) — and this function again satisfies periodic boundary conditions, 
so that -ffc^o (t) (sinh a) = — s^sinha. This function has two nodal regions, hence 
— must equal either the second or third eigenvalue of Hoa(t). On the other hand, 
we also know that L r d s a = 0, and this function is again periodic, so is also in the 
spectrum of Hq^(t). Since d s cr again has only two nodal regions, must also be 
either the second or the third eigenvalue. Putting these statements together shows 
that — is the second eigenvalue and is the third eigenvalue of Ho } q(t). 

To obtain the bounds on the bottom eigenvalue, note that the potential in L T .o 
satisfies the estimates 

-s 2 t (t 2 + 2) < -s 2 t t 2 cosh(2cr) < -s 2 t t 2 

since 

t 2 cosh(2cr) = r 2 + 2 t 2 sinh 2 a, 

and t 2 sinh 2 a < 1. The estimate for Ao(t, a) is then straightforward by mono- 
tonicity, but cf. below for the precise form of the eigendata of these a-quasiperiodic 
problems when the potential is replaced by a constant. □ 
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We have now shown that 

BoWUBi(r) = [Ao(r ) 0),Ao(r,7r)]U[Ai(r ) 7r),-4] C (-00, 0), 

and 

Ba(r) = [0,A 2 (T,7r)]. 

Recalling our earlier remark that no band reduces to a point, we have A2(t, 7t) > 
and hence Bk(r) C K + when k > 3. 

Remark 4.5. The main conclusion of this discussion is that only the bands 
Bq(t) and -Bi(r) lie in the negative half-line; all other bands are wholly contained 
in the positive half-line. 

It will also be necessary later to have more refined information about the be- 
haviour of the a-quasiperiodic eigenvalues as r — > — 00. The key observation is that 
L Ti q converges uniformly in this limit to L := — <? 2 — 1, and hence the spectrum 
of L T fl converges to that of L. Clearly spec (L) = [— l,oo), but although L has 
constant coefficients we may still perform the Bloch wave analysis to decompose 
this spectral ray into an infinite union of spectral bands. In fact, the solutions of 
Lu = \u with u(2n) = e M u(0), u'(2vr) = e ic V(0) are given by e i « Q / 27r ) ±fc ) t , and so 
for a £ [0, 7r], we have 



A fc (o) 



2 



^- h k I — 1 when /.' is even 



(£-*) 



2 



1 when A: is odd. 



and the corresponding eigenfunctions are given by 

5 l (37-fc)* when k 



is even 

e i(^+k)t w hen k is odd. 



This gives the band structure 

B Q = [-1,-3/4], B x = [-3/4,0], B 2 = [0,5/4], etc. 

The apparent lack of smoothness in the band functions Afe(a) at a = 0, tt is due to 
the absence of gaps between these bands. 

We now use Lemma 3.2 to perturb off this limiting situation. 

PROPOSITION 4.6. For a e [0,tt], let (<p k (a), Afc(a)) and ((pk(r, a), Afc(r, a)) 
denote the eigenfunctions and eigenvalues for L T and L T 0; respectively. Then for 
t sufficiently negative, 

(4.5) A fc (r,a) = A fc (a)+0(|T|- 2 ), 9 T A fc (r, a) = 0(|r|- 3 ), 
and moreover, 

(4.6) ^ fc (r,a) =^ fe (a)+0(|r|- 2 ) 

uniformly along with all derivatives on < t < 2ir. 

Proof. The fact that the limits of Ajt(r, a) and ipk(T,a) converge to A^(a) 
and fk{oL) is clear from general theory. To get the estimate on the eigenvalue, use 
Lemma 3.2 to get that 

, r s 2 T 2 cosh(2a(s T i)) = 1 + 0(|t|- 2 ), and 

( > d T (s 2 T T 2 cosh(2a(s T t))) = 0(|r|^ 3 ). 
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Assuming that the functions ipk (t, a) are normalized to have L 2 norm equal to 1 , 
we use the standard formula from eigenvalue perturbation theory 

d T X k (r, a) = (d T L Tfi <p k (r, a), <p k (r, a)) L2 = C(|t|" 3 ), 

whic h follows directly from (4/7). Integrating from — oo t o r gives the first part of 
(4.5); then a standard perturbation argument yields (4.6). □ 

4.4.3. Spectral flow of H a (r). We now let H a (r) denote the operator L T acting 
on a-quasiperiodic functions. (Thus it is rotational invariant, and its reductions 
to the eigenspaces of the cross-section are the operators Hj, a {r).) The motivation 
for this whole paper is the fundamental observation that H a (r) has a nontrivial 
spectral flow as r decreases from to — oo. To see this, recall that we have already 
determined that its eigenvalues are of the form 

\ k (T,a) + s 2 T j 2 , j,k = 0,1,2,... 

In addition, when a £E [0, 7r], Ao(t, a) stays strictly negative and uniformly bounded 
away from zero as r — > — oo, and the same is true for Ai(r, a) for a S (0, tt] (but 
note that Ai(r, 0) = — s 2 —►()). In particular, we see that Ao(r, a) + s 2 j 2 is positive 
for all j (because of (3.5)), but eventually becomes negative for each fixed j. Thus 
more and more eigenvalues, corresponding to higher and higher eigenmodes on the 
cross-section, change sign from positive to negative. 

The following two propositions make this more precise. The first is valid for any 
j > 2 and for 'intermediate' values of r. The second gives a much more accurate 
estimate, but is only valid when \t\ is large. 

Proposition 4.7. Let j > 2. If the interval B (t) + j 2 s 2 contains 0, then 
Bq(t) +4j 2 s 2 is entirely contained in the positive axis (0,oo). 



Proof. It follows from Proposition [4.4| that 

Ao(r,a) > -(t 2 + 2)s 2 t . 

Hence, Bo(t) + j 2 s 2 cannot contain if r 2 < 2 and j > 2. Therefore, we can 
assume that r < — 

Now, suppose the result were to fail. Then for some j we would have both 

A (t,7t) +j 2 s 2 T > and A (r, 0) + 4j 2 s 2 < 0. 

This gives 

A (t,0)-4Ao(t,7t) <0. 
On the other hand, from Proposition 



-(r 2 + 2)s 2 T < A (r, 0) < A (r, tt) < | - r 2 s 2 , 

which implies that 

A (r, 0) - 4A (r, tt) > -1 + 3r 2 s 2 - 2s 2 . 
Combining these two inequalities we get 

si < 



3r 2 



But from (3.4), 



sl> 



1+T 2 
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These two inequalities are incompatible when r < — y2. d 

4.5. The index of L T . We now apply the information we have about the 
spectrum of H a (r) to the spectral flow and index of L T restricted to to the subspace 
functions which are T JjQ -symmetric and also invariant under the reflection (t, 6) t— > 
(— t, —6). Equivalently, we let L T act on functions u which satisfy 

u(t + 2n,6) =u(t,6 + a), u(t, 6 + 2n/j) = u(t, 6) and u(-t, -6) = u{t, 0) 

for all (t, 6) G [— 7r,7r] x S 1 . As before, we denote the H 2 closure of this space by 

Observe that the Fourier decomposition of any function u G H 2 a reads 
u(t,8)= J2 u n(t)e" 1 ' 6 

nSZ-{0} 

where u n satisfy u n (t + 2ir) = e" J " u n (t) and u n (— t) = — M_„(i). 

It is easy to see that the spectrum of L T acting on the space above is given by 



spec(L T )= |J |J (s 2 T n 2 j 2 + B k (T)) 



k>0 n>l 

Now, for j > 2, t S V~ and |a| < 7r/j, define /j jQ ,(r) to be the number of 
negative eigenvalues of the operator L T acting on H 2 a . We have established that 
there exists ar, < such that Ij^ a (r) = when r G (t*,0). As t decreases, is 
the first value at which an eigenvalue of H(t) crosses zero; since we are restricting 
to some proper subspace of functions, the spectrum of i? J;Q (r) may remain positive 
on a larger interval. Thus define 

Tj. a = inf {r : Ij, Q (r') = : r' G (r,0)}; 

hence all eigenvalues of Hj^ a (r) remain positive when r G (rj. Q ,0). So, for all 
t G (Tj iQ .), the index of L T is and for some r which is slightly smaller than 
Tj ta , the index of L T is at least 1. Observe that Tj l0l is a root of the equation 
\o(T,ja) + s 2 j 2 = 0. The next result shows that we have a good control on Tj jCt 
when j is large enough. 



Proposition 4.8. There exists a j > 2 such that if j > jo, and if \a\ < n, 
then 



T jt * = - / ^_ = + o(r 1 ). 



y/l - H2tt) 2 

Moreover Ij^ a (r) — for r > Tj >a and Ij t0t (r) > 1 for all r < Tj, a . 

Proof. The number Tj >a corresponds to the value of r for which Aq(t, j a) + 
s 2 j 2 = 0. The estimate here follows from (4.5) and the explicit expression for 



Afc(a). The point of requiring j to be large is that this forces r also to be large, 
and then we can use these asymptotics results. □ 

For general j > 2 our control on Tj <a is weaker. The proposition below is 



immediate from Proposition 4.4 



Proposition 4.9. For j > 2, we have Tj-, ,Tj G (—j,—\/j 2 — 2). 



The next result is a consequence of the proof of Lemma 3.2, and is the basis 



for proving that when j is large enough, the eigenvalue crossings are transversal. 
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Proposition 4.10. Fix a € [— ir, ir]. There exists j > 2 such that if j > jo 
and if Afc(r, a) + s 2 j 2 = 0, k = 0, 1, then 

d T (X k (r,a) + s 2 T j 2 ) < 0. 

PROOF. First, from Lemma |3l^, we know that d T Xk{r, a) — 0(|t|~ 3 ). Also, 
by f!|), s 2 T = 1/r 2 + C(|r|- 4 ), and d T s 2 T = -2/r 3 + 0(|t|^ 5 ). Hence 

d T (X k (r,a) + s 2 T f) = -2j 2 r- 3 + 0(\r\- 3 ). 

Then observe that, since Afc(r, a) + s 2 j 2 = (recall that is bounded since we 
only need consider k = 0,1), then s r ~ j -1 ~ |t| _1 , so 9 T (Afc(r, a) + s^j 2 ) is 
certainly negative when j is large enough. □ 

5. CMC deformations of nodoids 

We now employ the preceding results about the Jacobi operator to deduce the 
existence of CMC surfaces with Tj-. Q symmetry. 

5.1. The mean curvature operator. For any t£P, parametrize D T by 

X T (t,6) := - (re^O^ cos 0, r e <M s -*) sin0,K T (M)) . 

The unit normal at X T (t, 0) is then 

(5.1) N T (t,9) :— (t smha T (s T t) cos0,r sinher T (s r i) sin 9, — d s a T (s T t)) ■ 
If w is any function which is C 2 small, then let D w the image of the map 

X w = X T + wN T . 

Note that if w is Rj or S a or T JiCt symmetric, then D w has the same symmetries. 

A rather complicated nonlinear elliptic equation determines when D w has mean 
curvature 1. We write it in abbreviated form as 

(5.2) L T w + Q T {w) = 0, 

where L T is the multiple of the Jacobi operator which we have been studying, and 
Q T is a second order nonlinear differential operator which vanishes quadratically at 



w = 0. More precise information about the structure of Q T is given in |11|. Note 
that in the (t, 9) coordinate system, Q T is 2 n periodic. 

Because this is a nonlinear problem, we shall use the function spaces Cj'£ (M x 
S 1 ), defined for any j > 2 and a € [—n/j^/j] to contain all C k ^ functions which 
are Tj jCe -symmetric and invariant under the reflection (t, 9) > (—t, —9): 

(K x S 1 ) = {u S C k ' (R x S 1 ) : u(t + 2ir, 9) = u(t, 9 + a), 
u{t,9 + 2n/j) = u{t,6) and u(-t, -6) = u(t, 6)} . 

Clearly 

L T :Cf;f(Rx S l ) — .C°f(Rx S 1 ) 

and 

Q T : C|;f(M xS 1 )^ C°f (R x S 1 ) 

are smooth. 
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Consider the quotient of ]R 3 by the screw motion S a along the z-axis, with 
translation length 2ir; this is a 2-plane bundle over S 1 with holonomy a. D T /S a is 
a compact submanifold of this space, and we shall construct the surfaces bifurcating 
from D T as perturbations of D T /S a . 

5.2. Bifurcations. We are now in a position to prove the existence of families 
of (immersed) CMC surfaces which bifurcate off the Delaunay surfaces D T . 



The proof of Theorem 



1.2 follows from the general bifurcation theorem of 



Smoller and Wasserman [17], cf. also Theorem 13.10 in |16|. To apply this re- 



sult, we require only the fact that, by definition, the index Ij ta (r) is for any 



t > Tj t0l and, by Proposition 4.7, the index Ij, a (r) is exactly 1 for some r < Tj >a , 
but close to Tj <a . This change of multiplicity ensures the existence of a bifurcation 
for the nonlinear problem (5.2) in the space C?'f(Rx5 1 ) (modulo the screw motion 
S a ). When a = 0, this produces nonrotationally invariant CMC surfaces which are 
periodic and Rj -symmetric. When a^O this produces nonrotationally invariant 



CMC surfaces which are T Jia -symmetric. By Proposition 4.9 we get the estimate 
on the location of this bifurcation point. 

One defect of this general theorem is that we obtain no information about 
whether this bifurcation gives a smooth connected branch of solutions. For this wc 
require nondegenerate crossing of a simple eigenvalue. However, in Proposition 4. 10| , 



we have verified this hypothesis when j and hence r is sufficiently large. Therefore, 
we obtain Theorem 1.3 from the theorem of Crandall and Rabinowitz [E|, cf. also 



Theorem 13.5 in [16 



We conclude by noting that there are many unresolved questions concerning 
the surfaces we have produced. The most obvious one concerns the existence of 
a second bifurcation which arises when Ai(t, a) + s^j 2 crosses zero. This should 
not be difficult to obtain, and requires only a slight elaboration of the techniques 
and estimates we have been using. One complication here is how to separate the 
second bifurcation for some smaller value of j occurring at the same value of r as 
the first bifurcation for a larger value of j. We have not pursued this because our 
information about these surfaces is so limited, since we have proved their existence 
using an abstract functional analytic technique. The most interesting problem is 
to globalize this construction and find a complete characterization of all immersed 
cylindrically bounded CMC surfaces with two ends. 
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